Liquid ground state, gap and excited states of a strongly correlated spin chain 
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We present an exact solution of an experimentally realizable and strongly interacting one- 
dimensional spin system which is a limiting case of a quantum Ising model with long range in- 
teraction in a transverse and longitudinal field. Pronounced quantum fluctuations lead to a strongly 
correlated liquid ground state. For open boundary conditions the ground state manifold consists 
of four degenerate sectors whose quantum numbers are determined by the orientation of the edge 
spins. Explicit expressions for the entanglement properties, the excitation gap as well as the exact 
wave functions for a couple of excited states are analytically derived and discussed. 



PACS numbers: 67.85.-d,75. lO.Kt, 75.10. Jm,05.30.Rt 

In low dimensional systems strong quantum fluctu- 
ations can inhibit the formation of long range order. 
A particularly fascinating class where this is the case 
are spin liquids [1]. Very recent numerical theoretical 
work has revealed and explored spin liquid phases in 
two-dimensional systems with possible experimental re- 
alizability - the anti-ferromagnetic Heisenberg model on 
a Kagome lattice [2] and the frustrated XY-model on 
the honeycomb lattice J3^. An extensively studied one- 
dimensional system which exhibits spin liquid behavior 
and which is amenable to analytic treatment is the cel- 
ebrated spin 1 chain due to Affeck, Kennedy, Lieb and 
Tasaki (AKLT) This model has no free parameters 
and the wave function of the ground state, which shows 
short-ranged entanglement and a hidden string-order, is 
known analytically. Initially studied in the context of 
valence bond solids the ground state has been demon- 
strated to be also of practical relevance e.g. as a resource 
for measurement based quantum computation [5^ . More- 
over, the AKLT model has provided valuable guidance for 
spin models that cannot be treated analytically but are 
located in some proximity to it in the parameter space, 
such as the spin 1 Heisenberg chain [5HH] or certain spin 
1/2 ladder systems [Hllin]- In spite of this great success 
and the detailed knowledge of the AKLT wave function 
[TTl W]\ the analytical construction of excited state wave 
functions or an exact calculation of the energy gap has 
remained elusive. 

In this work we explore a one-parameter spin 1/2 chain 
that is a limiting case of a quantum Ising model. Like the 
AKLT spin model it belongs to the class of frustration- 
free Hamiltonians 13 and has a quantum liquid ground 
state with a simple matrix product state (MPS) repre- 
sentation [13] . We present a detailed analysis of its corre- 
lations, degeneracy and entanglement properties. More 
importantly, we provide analytical expressions for the en- 
ergy gap and discuss the construction of the exact wave 
functions of a few excited states. Having access to these 
quantities is uncommon in non-trivial frustration-free 
systems with non-commuting local Hamiltonians. Here 
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FIG. 1. (a) Representation of the ground state in terms of 
classical spin configuration (here for L = 6 spins) which can be 
translated into dimer arrangements of an interstitial lattice. 
In general all accessible configurations contribute. The weight 
of each configuration is determined by the parameter z. (b) 
Chain of L spins with open boundary conditions. Here the 
ground state is four-fold degenerate and the individual sectors 
are labeled by the orientation of the two edge spins (v, /j.) with 
Ujfi =ti (only three sectors are shown). Depending on the 
state of the edge spins the number of spins contributing to 
the bulk differs. 



the existence of an energy gap can be usually established 
[T5] but neither its exact value nor the structure of ex- 
cited states are known. We believe that our model is 
particularly appealing for it is experimentally realizable, 
possesses a ground state with a simple MPS wave func- 
tion and has a tuneable excitation gap. At the same time 
it is amenable to analytical treatment such that proper- 
ties that go beyond the ground state can be exactly de- 
rived. It can therefore serve as a starting point for the 
analysis of other strongly correlated systems that depart 
from the exactly solvable parameter space. 
Hamiltonian — The spin model we are considering here 
consists of a one-dimensional chain of L spin 1/2 particles 
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that interact via a three-body interaction. The Hamilto- 
nian is given by 



H = Hn 



with 



L-1 



Ho = J2 Pk-i H + zPk + z-^nk\ Pk+i 



k=2 



L-1 

k=2 



(1) 
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Here P^ = (1 — cr^)/2 is the projector on the down- 
state of the fc-th spin cr^ and cr^ are Paiih 
spin matrices and the number operator = f — 
Pfc is the complement of Pk- Local Hamiltonians hk 
that belong to adjacent sites do not commute. iJb 
contains the boundary terms and reads for periodic 
boundary conditions ffb = Pl [c^ + z~^ni + zPi] P2 + 
Pl-1 [tr^ + z~^nL + zPl\ Pi = hi +11^. In this case the 
Hamiltonian is symmetric under inversion and transla- 
tions of the lattice sites and depends solely on the pa- 
rameter z which we take to be positive and real. In Ref. 
[16j it was explicitly shown how this model can be physi- 
cally realized within a lattice gas of cold atoms but imple- 
mentations with polar molecules [171 118| or trapped ions 
[TO] are in principle equally possible. Those systems are 
governed by a Hamiltonian which maps on the quantum 
Ising model in a transverse and longitudinal field with 
nearest and next-nearest neighbor interaction [5D]: 

L 

Hph = ^ [(J^ + f{z) Uk + v UkUk+i + z nknk+2] ■ 
fe=i 

Hamiltonian ([ij is a special case that emerges in the 
limit w 3> 1 (exclusion of neighboring excitations) and 
with f{z) = z~^ — 3z. In the above-mentioned physi- 
cal realizations the parameter z can be changed through 
an adjustment of experimental parameters such as the 
strength of laser driving and the two-body interaction po- 
tential. Hamiltonian ([T]) possesses the L conserved quan- 
tities TifcTifc+i, i.e., [nkUk+i, H] = OVfc. In this work we 
are interested in the subspace in which all of these opera- 
tors have eigenvalue zero (zero-subspace) , i.e. an up-spin 
is always accompanied by a down-spin on either side. The 
previously mentioned physical implementations of H are 
naturally confined to this subspace whose dimension is 
(f>^ with (f> = {1 + \/5)/2 being the golden ratio. 
Ground state — The ground state of Hamiltonian ([T]) 
obeys h^ \z) = QVk, i.e., it is the ground state for each 
of the positive semi-dcfinitc local Hamiltonians /ij. in eq. 
([1]). This is the defining property of a frustration- free 
system [T^]. The ground state energy is hence zero and 
the corresponding wave function reads explicitly 



l[Ali-z)\U...i). 



(3) 



Here the operator AI,{z) = exp [zPk-ia^Pk+i] with 
cr+ = {a^ + i<Jy)/2 creates a spin in the state jj,) -f z |t) 



on the fc-th site when applied to the fully polarized state 
III ... I). The projection operators in the exponential of 
A^(z) ensure that no two adjacent spins are simultane- 
ously in the up-state and one thus remains in the zero 
subspace with respect to the operators rikUk+i- Further- 
more, one can easily verify that [Aj, (x) , A\^ (y)] = OV k,m 
and Al{x)Al{y)^Alix + y). 

The state ([s]) can be written as a weighted superposi- 
tion of all possible dimer arrangements where an up-spin 
is identified as a dimer occupying an interstitial lattice 
as shown in Fig. [iji. The ground state is thus a su- 
perposition of a huge number of classical spin configura- 
tions as is indicated in the Figure. In the extreme case 
z = 1 all spin configurations contribute with equal weight. 
This is a manifestation of the strong quantum fluctua- 
tions that probe the entire accessible Hilbert space and 
prevent order from being formed. In the dimer picture it 
becomes evident that the normalization constant of the 
state ([3| is equivalent to the partition function of a one- 
dimensional gas of hard dimers at fugacity pjjj, and 
hence N{z) = [(1 + Vl + 4z2)/2]^. 

Let us continue by discussing some properties of 
the state It has a MPS representation 2lJ |z) — 

[N{z)]-'/^J:^,,..^,=l.t^^[XnX,....X,A\H.^2,'■.■,^L) 
with the matrices X-^ = and X_i = P — za^ . The 
bond dimension of this MPS is two and entanglement 
is only present between nearest neighbors [22]. The 
two spin reduced density matrix Pkk+i in the basis 
{lt>fc i;>fe+i , i;)fe lt)fc+i , i;)fe \l)k+i} is given by 



(n) C{z)/2 -{n)/z 
Pkk+i = I C{z)/2 {n) - {n) /z 
- (n) /z -{n) /z 1-2 (n) 



(4) 



Here we have abbreviated the density (n) = 
[1 — l/Vl + 4z2J /2 and the concurrence C{z) — 
2 (n) [2 (n) - 1] / [(n) - 1]. The concurrence < C(z) < 
1 is directly related to the entanglement of formation [23] 
but also is by itself a measure of entanglement, e.g. the 
fully entangled singlet state (l/-\/2)[|ti) — lit)] has a con- 
currence of one while an unentangled product state has 
zero concurrence. In the case of the density matrix Q 
C{z) assumes its maximum value C(zmax) = 6 — 4-\/2 = 
0.34 at z„iax = Vl + V2/V2 = 1.1. 

Spins that are separated by one or more sites are not 
entangled and are described by product states. However, 
one has to distinguish here between a separation by an 
odd/even number of sites. The corresponding density 
matrices in the standard basis read p°^'^ = [pf^ (8) p^- -I- 



> p„]/2 and p! 



ieven 



'P+]/2 with 



G, 



[{n) T G.,] /z 



- [(n) T G.,] /z 
1 - in) - G„- 



(5) 



Here Gij = \J\dij\ where dij = 

-2z^ 



[zV(l + 4z2)][(Vl+4z2. 
nected density-density correlation function. From these 



UiUj) - [Ui) {n.j) = 
l)/(2z2)]l'"j| is the con- 
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exponentially decaying correlations we can read off the 
correlation length ^: 



log 



1 - VlT4~2 

2^2 



(6) 



This shows that only for z ^ oo long range order (with- 
out long range entanglement) is present. 
Boundary conditions — So far all results were given for 
periodic boundary conditions, i.e., in the presence of 
-ffb: where the ground state is unique. This changes 
for free boundaries in which case iJb = 0. Here the 
Hamiltonian commutes with the operators <jI and erf', 
i.e. the magnetization of the edge spins is conserved. 
This leads to four disjoint sectors which can be la- 
beled by the magnetization of the edge spins {v, fi) 
with /i =tJ, as shown in Fig. [T|d. In the sec- 
tor (J,, J,) both edge spins are in the state |^) and we 
can replace the term h2 — Pi [a^ + zP2 + 2^^112] P3 
in Hamiltonian ^ by [cr^ -I- ZP2 + z^^n2] P3 and sim- 
ilarly for /iL-i- The ground state of the spin chain 
with fixed edges is again unique, has zero energy and 
can be written in the form ([3| with the difference that 
A\{-z), A\{-z) 1 and Ali-z) exp[-zCT^P3], 
AI^_^{—z) — > exp [— zPl-2ct:^^^] . Hence the ground 
state in this sector is = ||)^ |bulk(z))^^ ig) 14,)^ and 
its normalization constant is determined by the bulk wave 
function Nii{z) = ^^(bulk(z) | bulk(z))^^. This constant 
is equivalent to the partition function of hard dimers 
placed on a line with L — 2 sites and fugacity . Similar 
arguments lead to the construction of the ground state 
wave functions in the remaining three sectors. However, 
due to nin2 — ulul-i = a spin-up state at the edge 
signifies that the adjacent spin has to be in the down- 
state which reduces the number of spins contributing to 
the bulk (Fig. The fact the four degenerate ground 
states are gapped (discussed below) and that they can be 
distinguished by the edge spins might make them useful 
for the storage of quantum information. 
Spectrum — Let us return to periodic boundary condi- 
tions. We saw in eq. ([6| that the ground state exhibits 
density-density correlations that in general decay expo- 
nentially with the distance. Only in the limit z — 00 
the correlation length reaches the entire system length. 
This behavior indicates the existence of an excitation gap 
which closes as z approaches infinity. That this is indeed 
the case is shown in Fig. [2] in which displays numerical 
data for the energies of the ground state and the first 
nine excited states of a chain with 16 spins as a func- 
tion of z. The ordered ground state at z — >■ 00 is given 
by the symmetric superposition of the degenerate (anti- 
ferromagnetic) states \u) = Ititi •■•) and \d) = |itit •••)■ 
It is now important to understand whether the gap at 
finite z persists also in the thermodynamic limit. To get 
a first answer we have numerically calculated the spec- 
trum at 2: = 1 and varied the system size. Fig. [3] shows 




FIG. 2. Ground state and the first nine excited states (some 
of which are degenerate) of a chain of L = 16 spins. The gap 
between the ground state (blue) and first excited state (red) 
is finite except for z — > 00. The position of the red squares 
show the energy of the first excited state according to the 
analytical formula (|9|. 




FIG. 3. Ground state energy and the energy of the first three 
excited states as a function of the system size L. The data is 
calculated for z = 1 and periodic boundary conditions. If L is 
even the gap is system size independent and remains constant 
at a value of 2 — \/2 as predicted by eq. 



the corresponding data. The gap persists for all values 
L shown in the Figure while the separation between ex- 
cited states decreases. This suggests that excitations are 
gapped for all L and that the spectrum above the gap 
becomes continuous as L — > 00. Another observation is 
that for even values of L the excitation gap is independent 
of the system size. This feature is the decisive hint for 
the following analytical construction of the wave function 
of the first excited state. 

Excited states — We seek to construct excited states 
of the form \E) = X \ z) where X is an operator that 
creates an excitation on the ground state. We choose the 
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following ansatz for this operator: 



and 



X = 



(-1)" [arijn + P nm-irim+i] 



(7) 



where a and /3 are real numbers. This form is motivated 
by the following two observations: 

(i) The action of off-diagonal operators on the ground 
state, is proportional to the action of diagonal operators: 
By construction we know that hk\z) — which is the 
property of a frustration free Hamiltonian. Hence also 
nkhk \z) — and Pkh^ \z) — 0. Writing out the hk ex- 
plicitly and utilizing that ni^rik+i \z) — one finds the 
following relations 

Pk-ialPk+i \z) = -z~^nk \z) 
Pk-icrtPk+i \z) = -z[l- Hk-i - Uk 

-rik+i + nk-iUk+i] \z) 

This shows that an ansatz for X which contains number 
operators and products of number operators is already 
the most general one. 

(ii) The first excited state must be antisymmetric under 
a translation by one lattice site: In the limit 2 — >■ cxd 
the state ([3) reduces to the symmetric superposition 
\z — >■ 00) = [\u) + \d)] /V^, but also all spin configura- 
tions without three adjacent down-spins, i.e. |... Hi ...), 
are equally ground states of the system as their energy 
tends to zero (see Fig. [2]). The question is now which 
superposition of these configurations becomes the first 
excited state as z assumes a finite value. To see this, we 
perform perturbation theory at large but finite z. Here, 
the term of the Hamiltonian ([l]) that contains consti- 
tutes a perturbation that gives rise to virtual spin-flips, 
e.g. |...|ti---) \... Hi ■■■) |...|tl---)- A state of 
the ground state manifold (at z — > 00) with m up-spins 
can undergo m such virtual spin flips each of which low- 
ers the energy by — That means the more up-spins a 
configuration contains, the lower its energy will be. Since 
this state also has to be orthogonal to |z — 00) the an- 
tisymmetric superposition [\u) — \d)] / \/2 appears to be 
the correct first excited state in the 2: — )■ 00 limit. The 
antisymmetric behavior under shifts by one lattice site 
must be retained for all z as states of different z are adi- 
abatically connected. 

Using now the ansatz Q and calculating explicitly 
H\E) = HX\z) = E\E) = EX\z) leads to a closed 
set of equations for the coefficients a and (3 when L is 
even: 



- z 2z 
-z 2z + 2z-'^ 



E 



The corresponding eigenvalues and eigenvectors are 



-E4 



3 + z2 ± Vl + 6z2 + 
2z 



(8) 



(9) 



ip Vl + 6z2 
2z2 



+ 2^ 



,(10) 



respectively. E- corresponds to the energy of the first 
excited state and therefore to the excitation gap. The 
perfect agreement of this analytical expression with the 
numerical data is shown in Fig. [2j 

Let us finally analyze the structure of the first excited 
state more closely. We can express its (unnormalized) 
wave function in terms of the operators [x) : 

L 

\E-) = 5](-ir<_i(7-)AL+i(7-)k) (11) 

m— 1 

with 7_ = (2z/3)/(q;_). This form suggests a pair of ex- 
citations that travel on the ground state with a lattice 
momentum tt. To get a further idea about the appear- 
ance of this state it is also instructive to inspect the spin 
configurations that participate to \E^). We remember 
that for the liquid ground state these were in general all 
accessible configurations (see Fig. [l^) . In contrast to this 
\E^) contains only configurations that differ in the num- 
ber of up-spins on the two sublattices, i.e. the sublattices 
formed by the sites with even/odd label. Moreover, due 
to the second term in eq. ([7|, configurations where the 
two sublattices differ also in the number of pairs of adja- 
cent up-spins are weighted different to those where this 
is not the case. 

Outlook — In this work we have analytically solved cer- 
tain aspects of the model ([!]), but further questions re- 
main, e.g. 'What is the nature of the (quasi-particle like) 
excitations?' or 'Can the model actually be solved en- 
tirely?' Concerning the latter, one can attempt to con- 
struct higher excited states by multiple applications of 
the operator X to the ground state. However, this ap- 
proach fails. It would be also interesting to see whether 
the treatment similar to the one presented here can be 
successful for other frustration-free models and/or in 
higher dimensions [241. 
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